Diffusion on a solid surface: Anomalous is normal 
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We present a numerical study of classical particles diffusing on a solid surface. The particles' 
motion is modeled by an underdamped hangevin equation with ordinary thermal noise. The particle- 
surface interaction is described by a periodic or a random two dimensional potential. The model 
leads to a rich variety of different transport regimes, some of which correspond to anomalous diffusion 
such as has recently been observed in experiments and Monte Carlo simulations. We show that this 
anomalous behavior is controlled by the friction coefficient, and stress that it emerges naturally in 
a system described by ordinary canonical Maxwell-Boltzmann statistics. 



I 

C 

O 
o 



> 

00 

in 
o 

m 
O 



-a 
c 

o 
o 



X 



PACS: 05.40-a, 68.43.Jk, 68.35.Fx 

Diffusion of atoms, molecules, and clusters on solid 
surfaces occurs in a number of modern technologies in- 
volving self-assembled molecular film growth, catalysis, 
and surface-bound nanostructures [lj. The study of the 
motion of small and large organic molecules 0, y], and 
of adsorbed metal clusters composed of tens and even 
hundreds of atoms 0, 0| , has led to the unexpected ob- 
servation that, as with single atoms [6(, long jumps may 
play a dominant role in these motions. 

Theoretical, numerical, and phenomenological discus- 
sions of surface diffusion have led to the clear understand- 
ing that jumps beyond nearest neighbors are ubiquitous 
in some parameter regimes 0. However, these studies 
focus on the fact that the motion is necessarily diffusive 
(which is the case at very long time scales). The possi- 
bility that jumps can be so long as to lead to superdiffu- 
sive motion over appreciable intermediate time scales is 
recognized as an interesting problem, but one in which 
Levy walks or flights Q are invoked as a model input. 
AlthoughLevy-walk-like behavior is clearly observed in 
Hamiltonian systems 0, llOf and in microcanonical sim- 
ulations [l(J, uM i ^ is generally believed that such a fine 
signature of chaos is fully smeared away by thermal fluc- 
tuations. In the present work we show that Levy- like 
statistics appear quite naturally over long time scales 
within the usual Langevin framework for underdamped 
motion in a periodic or a random potential. 

While a detailed analysis of surface diffusion requires 
extensive calculations (e.g., ab initio, or molecular dy- 
namics), even the most powerful currently available com- 
puters can not carry such calculations to anywhere near 
experimentally relevant time scales 01 • Moreover, cur- 
rent experimental probes of the topography of surfaces, 
scanning tunneling microscopy and atomic force mi- 
croscopy, are usually carried out at relatively high tem- 
peratures, which leads to additional difficulties for first- 



principles calculations. Therefore, simpler approaches 
are essential and valuable 0. 

We consider a generic model of classical particles mov- 
ing in a two-dimensional potential, under the action of 
thermal fluctuations and dissipation, the important con- 
trol parameter being the friction coefficient. In spite of 
the simplicity of the model, we find that it is able to 
reproduce the entire range of experimentally and compu- 
tationally observed phenomena, from superdiffusion all 
the way to subdiffusion. 

The equation of motion of a particle of mass m on the 
surface is 



-W(x/A) - /ix + £(t) 



(1) 



where A is the characteristic length scale of the poten- 
tial. The parameter fi is the coefficient of friction, and 
the (,i(t) are mutually uncorrelated white noises that 
obey the fluctuation-dissipation relation (£,i{t)£j(t')) — 
2/ikBTSij6(t — t'). We first consider the nonseparable pe- 
riodic potential 
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where Vo is the barrier height at the saddle points. 

Equation Q can be rewritten in scaled dimensionless 
variables, r x = x/X, r y — y/X, and s = y/Vo / mXt, leav- 
ing only two independent parameters, the scaled temper- 
ature T and the scaled dissipation 7, 



T = k B T/V , 7 = fiX/y/mVo. 



(3) 



We study four properties of the motion of the particle: 
the mean square displacement, the dependence of the 
diffusion coefficient on friction, the probability density 
function of displacements, and the velocity power spec- 
trum. 

Normal diffusive behavior is characterized by a lin- 
ear time dependence of the mean square displacement, 
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FIG. 1: Left: A trajectory for 7 = 1 over t = 20,000 time 
units. Right: A trajectory for 7 = 0.04 over t — 15, 000 time 
units. The period of the potential is A = 4. Note the different 
scales in the two panels. 



Non-diffusive behavior shows a different 





with a > 1 (< 1) for 



(r 2 (s)) ~ 

time dependence, (r 2 (s)) 
superdiffusive (subdiffusive) motion. In Fig. ^ we show 
typical trajectories obtained for two friction coefficients 
upon numerical simulation of the equations of motion 
with T = 0.2 (we use this value throughout). One (left 
panel) is for a large friction coefficient, and the particle 
follows typical diffusive motion characterized by short 
steps of length ~ A and frequent changes in direction. 
The other (right panel) corresponds to a small friction 
coefficient and clearly shows the preponderance of long 
(~3> A) tracks along cartesian coordinates. 

The evolution of (r 2 ) averaged over 1000 particles is 
shown in Fig. |3 for several friction coefficients. For very 
long times the motion is diffusive, as expected, but for 
small 7 and at intermediate times there is clear superdif- 
fusive ballistic (a — 2) behavior over several decades in 
time. This is reflective of the long straight stretches seen 
in the low-7 trajectory in Fig. ^ We stress that this 
behavior has emerged naturally and has not required ex- 
plicit insertion of any but ordinary thermal fluctuations 
in the model. 

Eventhough the motion of the particle may include 
long superdiffusive stretches, at long times the motion 
is necessarily diffusive. The dependence of the diffusion 
coefficient on the friction for small and for large 7 can 
be obtained analytically using the approximate relation 
D w (Z 2 )/2t, where (I 2 ) is the mean square size of a 
jump out of one well and into another, and r _1 is the 
mean jump rate (related to the familiar "mean escape 
rate"). In the overdamped regime, jumps typically occur 
from one well to a neighboring well, so (l 2 ) « 1. Familiar 
Kramers formulas can be used to obtain the mean escape 
rate 03, with the result 
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(4) 



The 7 1 dependence of D arises because (I 2 ) is inde- 
pendent of 7 while r ~ 7. In the underdamped limit, 
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FIG. 2: (r 2 )/4s for a particle in the periodic potential, for 
7 = 0.0004 (solid), 0.004 (dotted), 0.04 (dashed) and 0.4 (dot- 
dashed). The straight-line segment has unit slope as a guide 
to the eye. Inset: Diffusion coefficient as a function of 7. The 
solid lines correspond to the theoretical calculations (Eqs. (QJ- 
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FIG. 3: Log- log plot of probability distribution functions for 
the particle displacement for 7 = 0.0004 and three different 
values of time intervals : t s = 20 (squares), 100 (circles) and 
1000 (triangles), and for 7 = 1.0 at the time interval r s = 100 
(solid curve). 



standard results |l|| lead to r 
(I 2 ) ~ 7~ 2 Hj. One obtains 
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I/7 [14|, but now 



(5) 



that is, again an inverse dependence on friction. The 
theoretical diffusion coefficient as a function of the fric- 
tion parameter is shown as the solid curves in the inset 
of Fig. |21 The symbols are the simulation results. The 
7 _1 dependences have been noted in the literature 0, 
but we have provided explicit forms with no adjustable 
parameters. 

The diffusion coefficient characterizes only the very 



long time asymptotic dynamics. The interesting inter- 
mediate dynamics in the low friction regime that gives 
rise to long stretches of ballistic motion is reflected in 
the probability distribution function (pdf) of particle dis- 
placements r. This pdf is shown in Fig. for 7 = 0.0004 
and tree different time intervals t s . For comparison, we 
also show a typical pdf for high damping (7 = 1) at the 
intermediate time interval. In the high-7 curve the high- 
est maximum corresponds to no jumps (by far the most 
likely event at short times). The next is associated with 
jumps to a nearest neighbor well, and so on. In con- 
trast, the low-7 curves show a very different behavior, 
with features strongly resembling those of a Levy-walk 
model [flUl^: a peak at small displacements, a power-law 
intermediate regime, and a side hump at high displace- 
ments. Each of these is a distinct signature of Levy- walk- 
like dynamics, but one must be cautious in the detailed 
interpretation of these components. The persistent small 
displacement peak is associated with long trapping peri- 
ods during which a particle does not move at all because 
its energy is not sufficient to overcome the barrier. The 
high displacement peak, which moves outward with ve- 
locity of order unity, is associated with ballistic motion 
of those particles that acquire enough energy to move 
(and lose it very slowly). Genuine Levy- walk dynamics 
also exhibit a low displacement peak and a superdiffu- 
sive peak separated by a power law behavior, but there 
are some important differences. First, our distribution 
reflects ballistic transport in the intermediate regime (in 
the language of Ref. 15], ballistic transport occurs when 
< a < 1 in the Levy model) , whereas the regime where 
the Levy model shows the features we have described is 
associated with sub-ballistic (but still superdiffusive) be- 
havior (again, in the languag e of Ref. 0, the behavior 
when 1 < a < 2). Second, the slope in our power law 
regime (approximately 0.7) is not related to the expo- 
nent a in the mean square displacement as it is for the 
Levy walk (where the slope is 4 — a). Third, our side 
hump is strongly broadened whereas the side hump in 
the Levy- walk model is associated with motion at a single 
constant velocity. In our case the velocity varies accord- 
ing to the equilibrium Maxwell-Boltzmann distribution. 
Nevertheless, the qualitative features of our distribution 
track those of the Levy walk. Note that the existence 
of the pronounced side hump reflects the fact that the 
particles performing long steps ( "flights" ) are those with 
a velocity in the tail of the Maxwellian distribution. 

Long ballistic excursions imply velocity correlations 
over considerable time intervals. The velocity power 
spectrum S(uj) = (v(w)v(— u>)) for different values of 7 
is shown in Fig. 0] The pronounced peak at ujq = 
is associated with small oscillations in one well. At lower 
frequencies, u> <C u>o, one observes a power-law growth of 
S(ui) that corresponds to the persistent time correlations 
associated with the ballistic excursions. At even smaller 
frequencies there is a crossover to S(w) =const= D, in- 
dicating full decorrelation and pure diffusion. 

Disorder in surfaces occurs due to the presence of 
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FIG. 4: S(lo) for 7 = 0.0004,0.004 and 0.04 in decreasing 
order on the ordinate. 
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FIG. 5: Random potential of the same average potential 
height as in periodic case. 



vacancies and other defects. We have generated a 
random potential surface described by a Gaussian dis- 
tribution with a correlation function (y(x)V(x.')) = 
(e/27rA ,2 )exp(-|x-x'| 2 /2A' 2 ) (see [H M, MM f °r 
details). We set the intensity e = 100 and the character- 
istic length A' = 4 in our simulations. A typical surface 
generated with this algorithm whose average height Vb 
equals that of the periodic potential is shown in Fig. [3] 

The exponent a in (r 2 (s)) ~ s a at intermediate times 
shows the entire range of behaviors from subdiffusive to 
superdiffusive with changing friction. In Fig. we show 
(r 2 (s))/4s, averaged over 5000 particles, as a function of 
time, for several values of 7. In the overdamped re gim e 
we clearly observe subdiffusive behavior (a < 1) |l9j| . 
while superdiffusive (a > 1) behavior is seen for very 
small 7. The exponents a, calculated over the last decade 
of the time variation of the mean square displacement 
within our finite simulation times, are plotted in the inset 
of Fig. as a function of 7. Although the subdiffusive 
behavior is probably the true asymptotic behavior in the 
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FIG. 6: (r 2 )/4s for a particle in the random potential for 
7 = 0.0004 (solid), 0.008 (dotted), 0.04 (dashed) and 0.8 (dot- 
dashed). The straight-line segment has unit slope as a guide 
to the eye. Inset: Exponents a versus friction coefficient 7. 

overdamped case [l9j |. further theoretical and numerical 
support are needed to assess whether superdiffusion is 
the asymptotic behavior in the underdamped case. 

We summarize our findings. We have explored the be- 
havior of a particle in a two-dimensional potential de- 
scribed by ordinary Langevin dynamics under conditions 
of thermal equilibrium. In a periodic potential, in the 
underdamped regime, the motion of the particle includes 



a ballistic range that can extend over many decades of 
time. The pdf of the particle's displacements under these 
conditions shows a structure strongly resembling one for 
Levy walks. This may explain a number of observations 
involving superdiffusivc motion of organic molecules Q 
and atomic clusters ja] on surfaces without the need to 
invoke extraordinary fluctuations beyond the usual ther- 
mal description. The long-time behavior is diffusive in 
all cases, and we have been able to predict theoretically 
the dependence of the diffusion coefficient on friction over 
essentially the entire range of values of the friction pa- 
rameter with no adjustable parameters. The situation in 
a random potential is even more complex, and exhibits a 
wide range of subdiffusive to superdiffusive regimes. Fur- 
ther analysis of the random potential case, and a more 
extensive presentation of the periodic problem, will be 
detailed elsewhere pol |. 
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